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Abstract 

A lattice model with a spatial dispersion corresponding to a power-law type is sug¬ 
gested. This model serves as a microscopic model for elastic continuum with power-law 
non-locality. We prove that the continuous limit maps of the equations for the lattice 
with the power-law spatial dispersion into the continuum equations with fractional gen¬ 
eralizations of the Laplacian operators. The suggested continuum equations, which are 
obtained from the lattice model, are fractional generalizations of the integral and gradient 
elasticity models. These equations of fractional elasticity are solved for two special static 
cases: fractional integral elasticity and fractional gradient elasticity. 

PACS: 45.10.Hj; 61.50.Ah; 62.20.Dc 


1 Introduction 

The integration and differentiation of any arbitrary order have a long history since 1695 [1] 
- [4], The theory of equations with derivatives and integrals of non-integer orders [3j ;4j [5] 
is powerful tool to describe the behavior of media and systems with spatial non-locality and 
memory of a power-law type. Fractional differential equations have a vast applications in physics 
and mechanics and form an area called fractional dynamics m - ns]- The fractional calculus 
allows us to formulate a fractional generalization of non-local elasticity models in two forms: 
the fractional differential (gradient) elasticity models (weak power-law non-locality) and the 
fractional integral models (strong power-law non-locality). Fractional integral elasticity models 
are considered in |T6]J - [22] . 

Lattice models are very important in elasticity theory (see for example [23] -[28]). I* 1 this 
paper we suggest a lattice model with a power-law spatial dispersion as a microscopic model of 
elastic continuum with weak and strong power-law non-locality. For this purpose, we consider 
the models of lattices and their corresponding continuous limits by using the methods suggested 
in [29, [30] (see also |3T] [32, 33]). In [294 SQj we prove that the continuum equations with 
fractional Laplacian in the Riesz form mm can be directly derived from lattice models with 
different types of long-range interactions. In this paper, we show how the continuous limit for 
the lattice with a power-law spatial dispersion gives the corresponding continuum equation of 
the fractional integral and gradient elasticity. The continuum equations of fractional elasticity, 
which are derived from the lattice model, are solved for two static cases: fractional integral 
elasticity and fractional gradient elasticity. 
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2 Equations for Displacement of Lattice Particles 

Let us consider a lattice model where all particles are displaced in one direction. We assume 
that the displacement of particles from its equilibrium position is determined by a scalar held. 
It allows us to describe the properties of the lattice and its continuum limit by using simple 
equations. 

The equations of motion for a one-dimensional lattice system of interacting particles have 
the form 

M Q t 2 = 9 K i n i m ) (“n(f) - u m (t)j + F(n ), (1) 

m= — OO 

m^n 

where u n (t) = u(n,t ) is the displacement of n particle from its equilibrium position, g is the 
coupling constant for interparticle interactions in the lattice, the terms F(n) characterize an 
interaction of the particles with the external on-site force. For simplicity, we assume that all 
particles have the same mass M. The elements K (n, m) of equation (P) describe the interparticle 
interaction in the lattice. For an unbounded homogeneous lattice, due to its homogeneity, 
K(n,m) has the form K(n,m) = K(n — m). Equations of motion ([I]) is invariant with respect 
to its displacement of the lattice as a whole, provided external forces are absent. It should be 
noted that the noninvariant terms lead to divergences in the continuous limit [13] . 

3 Transform Operation for Lattice Models 

In order to define the operation that transforms the lattice equations for u n (t) into the contin¬ 
uum equation for a scalar held u(x,t), we use the methods suggested in |2D| 3D]: We consider 
u n (t) as Fourier series coefficients of some function u(k,t ) on [— k 0 /2, k 0 /2], then we use the 
continuous limit ko > oo to obtain u(k,t), and finally we apply the inverse Fourier integral 
transformation to obtain u(x,t). Diagrammatically this can be written in the following form. 

u n {t) uik,t) u(k,t ) ——» u(x,t) (2) 

We performed similar transformations for differential equations to map the lattice equation 
into an equation for the elastic continuum. We can represent these sets of operations in the 
form of the following diagrams. 


Equation for u n (t ) 


From Lattice to Continuum 
- > 


Equation for u(x,t ) 


Fourier series transform 


Inverse Fourier integral transform 


Equation for u(k,t ) 


-» 

Limit Ax —^0 


Equation for u{k , t) 


( 3 ) 


Therefore the transformation operation that maps our lattice model into a continuum model 
is a sequence of the following three actions (for details see [29] [30]): 

1. The Fourier series transform : u n (t) Ea{ u n {t)} = u(k,t ) that is defined by 


u{k, t ) 


Too 

Y, Mt)e- ikx ' 


n =—oo 


EA{u n (t)}, 


( 4 ) 
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i r+ko/2 

u n (t) = -j~ dk u(k,t) e lkXn = {u(k, t)}, (5) 

ko J~k 0 /2 

where x n = nAx, and Ax = 2ir/ko is the inter-particle distance. For simplicity we assume 
that all lattice particles have the same inter-particle distance Ax. 


2. The passage to the limit Ax —> 0 {k 0 ► cx)) denoted by Lim : u{k, t ) —* Lim{h(fc, t)} — 
u{k,t). The function u(k,t) can be derived from u[k,t) in the limit Ax —> 0. Note that 
u(k,t ) is a Fourier integral transform of the held u(x,t), and u{k,t) is a Fourier series 
transform of u n (t ), where we use 


/ \ 2tt . . 

Unit) = —U{x ni t) 
ko 

considering x n = nAx = 27m /ko —> x. 

3. The inverse Fourier integral transform : uik, t) —y £)} — n(x, t) is defined 

by 


u(k, t) = 


/* + OO 
' —OO 

/*+oo 


u(x ]t ) = - 


dx e lkx uix, t) = F{uix, t)}, 
dk e ikx uik,t) = F~ l {uik,t)}. 


( 6 ) 

(7) 


Using the suggested notations we can represent diagram (J2D in the following form, 


Unit) 


From Particle to Field / 
- > U(X,t) 





uik, t) 


-» uik A) . 

Lim Ax—>0 


( 8 ) 


The combination of these three actions J r_1 , Lim, and J-& allows us to realize the transfor¬ 
mation of lattice models into continuum models [29] [30]. 

Note that equations (j4j) and (J5]) in the limit Ax —> 0 (fc 0 —> oo) are used to obtain the 
Fourier integral transform equations (E]) and (171) , where the sum is changed by the integral. 

Let us give the statement that describes the Fourier series transform for the equations for 
displacement of lattice particles (jT]) . The Fourier series transform J -a maps the lattice equations 
of motion 

M = g ^ K(n,m) (u n (t) - u m (t)j + F(n), (9) 

m=— oo 
m^n 

where iL(n, m) satishes the conditions 


OO 

K(n, m) = Kin — m) = K (m — n), |/t(n)| 2 < oo, 

71=1 


into the continuum equation 

M d U Q^2^ = ^(-^(°) - KikAx)^ uik, t) + ^{Fin)}, 


( 10 ) 


( 11 ) 
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where u(k,t ) = K(kAx ) = J-A{i^(n)}, and is an operator notation for the 

Fourier series transform. Here we use the notation 

+oo 

K(kAx) = J2 e~ iknAx K(n) . (12) 

n=— oo 
0 

Using K(—n) = K(n), the function (fT21) can be represented by 

+OO +00 

K(kAx) = K ( n ) {e~ iknAx + e iknAx ) = 2 K(n) cos ( kAx ), (13) 

n=l n =1 

and 

K( 0 ) 

For details see [29, 30] and [ IB] . 

4 Weak Spatial Dispersion of Power-Law Type 

In lattice models, dispersion is associated with different properties of the wave such as its fre¬ 
quency, wavelength, wave-number, amplitude and others. Spatial dispersion is the dependence 
of the elastic waves on the wave vector. In the model that is described by equation (fl]) the 
spatial dispersion means the dependence of the kernel Ji(|k|) on the wave vector k. This de¬ 
pendence is caused by non-local interactions in the elastic continuum. The spatial dispersion 
leads to a non-local connection between the stress tensor Uki and the strain tensor Eki- The 
tensor <Jki at any point r of the continuum is not uniquely defined by the values of Eki at this 
point. It also depends on the values of Eki at neighboring points r', located near the point r. 
Qualitatively describing the process in a lattice with spatial dispersion implies that the fields 
of the elastic wave moves particles from their equilibrium positions at a given point r, which 
causes an additional shift of the particles in neighboring and more distant points r' in some 
neighborhood region. Therefore, the properties of the continuum, and hence the stress tensor 
held aki depend on the values of strain tensor held Eki not only in a selected point, but also in 
some neighborhood region. The size i? 0 of the area of the mutual influence are usually of the 
order of the interparticle distance in the lattice. The wavelength, A, of elastic waves is several 
orders larger than the size of this region, so for a region of size i? 0 the values of the held of the 
elasticity wave do not change. In other words, the wavelength A usually holds kRo ~ Ro/X 1. 
In such a lattice the spatial dispersion is weak. 

To describe the dynamics of the lattice it is enough to know the dependence of the function 
A'(|k|) for small values of k — |k| only. Therefore we can replace this function by the Taylor 
polynomial. The weak spatial dispersion in the media with a power-law type of non-locality 
cannot be described by the usual Taylor approximation. The fractional Taylor series can be 
very useful for approximating non-integer power-law functions |36j . This is due to the fact that 
the usual Taylor series for the power-law function has an infinite number of terms. Using the 
fractional Taylor’s formula we obtain a finite number of terms. 


+oo 




(14) 


71 = — OO 
0 


n =1 
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For an isotropic linear medium with the weak spatial dispersion the function K (|k|) can be 
represented in the form 


N 


A-(|k|)« jr(0) + £o < , j |kp, 

3 = 1 


(15) 


where the frequency dispersion is neglected, and thus K( 0), and a aj (j = 1 do not 

depend on the frequency uj. 

If a.j = j for all j G N, we can use the usual Taylor’s formula. In this case we have 
a well-known weak spatial dispersion. In general, we should use a fractional generalization 
of the Taylor’s series [3j [1J [37, [38], 3S SOI Slj- The orders of the fractional Taylor series 
approximation should be correlated with the orders of power-laws of weak spatial dispersions, 
which are experimentally determined. In this case the fractional Taylor series approximation 
of K (|k|) will be the best approximation. 

We consider models of lattices with weak spatial dispersion and their continuous limits by 
using the methods suggested in [29030]. In this limit we obtain the fractional Laplacian in the 
Riesz’s form since the inverse Fourier’s integral transform of |k| Q gives the fractional Laplacian 
(—A)"/ 2 . Note that some equivalence of lattice-type networks with long-range interactions and 
continuum models for elasticity theory is considered in [42], where the Marchaud fractional 
derivatives are used. 

For a > 0 and x G M n , the fractional Laplacian in the Riesz’s form is defined in terms of 
the Fourier transform T by 


((-A)«/ 2 /)(x) = /-‘(liftTTOM). (16) 

For a > 0, the fractional Laplacian in the Riesz’s form usually is defined in the form of the 
hyper-singular integral by 


((-a r ' 2 j )( X ) 


i 

d n (m, a) 


1 

2 j OL-\-n 


(A ?f){z)dz, 


where m > a, and (A ™f)(z) is a finite difference of order m of a function f(x) with a vector 
step z G M n and centered at the point x G 


Til | 

<A r/)w = D-i)‘«(^i/(*-**)■ 


The constant d n (m, a) is defined by 


_t T 1+n/2 A m (a) _ 

2“r(l + cx/2)r(n/2 + a/2) sin(7ro;/2) ’ 


where 


m 


A m (a) = > f-iy 1 - 


m\ 


3=0 


j\(m — j)\ 


Note that the hyper-singular integral ((—A ) a ^ 2 f)[x) does not depend on the choice of m > a. 
The Fourier transform T of the fractional Laplacian is given by 


(^(-a r> 2 /m = ifci“(w/)(fc). 
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This equation is valid for Lizorkin space [3111] and the space C'°°(M ri ) of infinitely differentiable 
functions on M n with compact support. 


5 Fractional Elasticity Equation from Lattice Model 


In the continuous limit the equations for lattices with weak spatial dispersion of a power-law 
type gives the continuum equation for the fractional elasticity model. 

Statement. In the continuous limit Ax —y 0, the lattice equations 

o2 / 1 \ +oo 

M U fo 2 = 9 ^2 K(n - m) (u n (t ) - u m (t)^ + F(n), (17) 

m=—oo 

m^n 

with weak spatial dispersion of the form [T5\) gives the fractional continuum equation 


d 2 u(x, t) 
dt 2 


U P 


(18) 


with the fractional Laplacian (— A) aj//2 of order otj. Here the variables x and Ax are dimen¬ 
sionless, /( x) = F{x)/A\Ax\, p = M/(A\Ax\), A is the cross-sectional area of the medium, 
and 


n ga a .\Ax\ c 

— - 


M 


(j = 1....,JV) 


(19) 


are finite parameters. 


Proof. Let us describe the main points in the proof of this statement. Here we use the 
methods from |29, 30] and ra- The Fourier series transform Fa of equation (ITTj) gives CD. 
After division by the cross-sectional area of the medium A and the inter-particle distance |Aa;|, 
the limit Ax —> 0 for equation dTTjl gives 

8 2 JA n\Ar\ a i 1 

W2 u(k, t) = J2 lL JJ L - W n(k, t) + -V A {/(n)}, (20) 

3 =1 ^ 

where p = M/A\Ax\ is the mass density, f(n) = F{n)/A\Ax\ is the force density, |Ax| is the 
inter-particle distance, A is the cross-sectional area of the material, and 

JCa jA (k) = -a aj \k\ a i. 

Here we use (lT5ll . with G aj (j = 1,...,N) are hnite parameters that are defined by (lT9j) . The 
expression /C q ,a(^) can be considered as a Fourier series transform of the interaction term. 
Note that ga a . —> oo for the limit Ax —>• 0, if G a . are hnite parameters. 

In the limit Ax —> 0, equation (120]1 gives 

AAA = G a , K a ,(k) u(k, t ) + (21) 

3=1 P 
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where lC aj (k) = Lim/C Q) . ; a(&) = — a aj \k\ aj . The inverse Fourier transform of (12(1 gives (fTSlh 
Here, we use the connection between the Riesz fractional Laplacian and its Fourier transform 
lasiE] in the form \k\ aj <—)■ (—A)" 3 '/ 2 . 


To illustrate this Statement we give a few examples. 

Example 1. If we can use the weak spatial dispersion in the form 

K(k) « K (0) + a 2 k 2 , 

then we obtain the well-known one-dimensional equation for elastic continuum 

9 ^ — = G 2 A U(X, t) + - f(x), 
ot~ p 


( 22 ) 


(23) 


where 


G 2 = 


ga 2 |Aaf 

MA 


E 

P 


Here E = K\Ax\/A is Young’s modulus, K = ga- 2 is the spring stiffness, and p = M/A\Ax\ is 
the mass density. The corresponding dispersion relation is uJ 2 {k) = G- 2 k 2 . 

Example 2. If the spatial dispersion law has the form 


K(k) « K( 0) + a 2 k 2 + a 4 A: 4 , 
then we derive the equation of the gradient elasticity 

d 2 u(x,t ) 


as 


dt 2 


= G 2 A u(x, t ) — G 4 A 2 u(x, f) -|— fix), 

P 


(24) 


(25) 


where the constant of phenomenology model is associated with the lattice constants by 

ga 4 |Aa ;| 4 a 4 -E'|Aa ;| 2 


Ga — 


M A 


a 2 p 


The correcponding dispersion relation is thus ui 2 (k) = G 2 k 2 + G^k 4 . Note that the parameter 
l 2 of the gradient elasticity is related with the coupling constants of the lattice by the equation 


l 2 = 


|CI 4 1 |Ax|" 

I ci 2 1 


(26) 


The sign of the second-gradient term is defined by sgn(a 4 /a 2 ). Note that earlier it was thought 
that a phenomenological model of the gradient elasticity with a minus sign does not have the 
appropriate microscopic model, and it is thus considered one of its main weaknesses [3T]. The 
proposed lattice model radically changes the situation. 

Example 3. If we can use the fractional spatial dispersion law in the form 


K(k) « A'(0) + a 2 k 2 + a a k a , 
then the fractional elasticity equation is 

- U i, 2 ’^ = G 2 A u(x, t ) - G a (-A ) a/2 u{x, t) + - fix), 
ot z p 


(27) 


(28) 
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where 


ga a \Ax\ a a a E \Ax\ a 2 

M A clo p 


a 2 p 


The correcponding dispersion relation is co 2 (k) = G 2 A ; 2 + G a \k\ a . Equation ()28l) define the 
fractional elasticity model for the one-dimensional case (1 6 1 ). 


6 Solution of Fractional Elasticity Equations 


Using the same methods as above, we can derive a general model of three-dimensional lattice 
with fractional weak spatial dispersion of the form 


N 


A-(k) = A-(0) + 5>« J |k|“J. 


(29) 


Then the continuum equation for fractional elasticity model has the form 



(30) 


where we use Cj as a new notation for the constants instead of G aj used in the one-dimensional 
case. Note that r and r = |r| are dimensionless. 

Static equation and its solution 

Let us consider the statics ( du(r,t)/dt = 0, i.e. u(r,t) = u{ r)) in the suggested fractional 
elasticity model. Then equation (13(1 gives 



Equation (1301) has a particular solution w(r) for the case ajy > 1 and 7 ^ 0 (see Section 


5.5.1. pages 341-344 in [5]). The particular solution is represented in the form of the convolution 
of the functions G”(r) and /( r) as follow 



(32) 


where n = 1, 2, 3 and the function G”(r) is the Green function that is given by 



where a = (aq,..., a^). The Green function (1331) can be simplified by using the relation (Lemma 
25.1 of [3j S]) °f the form 



e *( k ' r ) /(|k|) d n k 



(34) 
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Here J v is the Bessel function of the first kind. As a result, the Fourier transform of a radial 
function is also a radial function. Using relation (1341) . the Green function (|33|) can be represented 
(see Theorem 5.22 in [5]) in the form of an integral with respect to one parameter A, 



where n = 1,2,3 and a = (cq,..., a m ), and J( n _ 2)/2 is the Bessel function of the first kind, 
which can be represented as J±/ 2 (z) = \j2j-KZ sin(^) for the 3-dimensional case. 


Thomson’s problem for fractional integral and gradient elasticity 

If we have the dispersion law in the form 



JL( k ) m AT(0) + ajk ® + a 2 k 2 , 

(36) 

where 

a > 0, then we obtain the fractional elasticity equation 



c 2 An(r) - c a ((-A)" / 2 w)(r) + - /(r) = 0, 

P 

(37) 

where 

E ga 2 Ax 2 ga a \Ax\ a 

C2 ~ p~ M ’ Cq “ M ' 

(38) 

If a = 

4 we have the well-known static equation of the gradient elasticity |34j: 



c 2 Att(r) - c 4 A 2 -u(r) + - /(r) = 0, 

P 

(39) 

where 

_ l ,2 e ga±\Ax \ 4 

Ca — it t — 

p M 

(40) 


The second-gradient term is preceded by the sign that is defined by sgn(ga 4 ), where g a 2 > 0. 
Equation ()37l) with n = 3 has the particular solution j5j of the form 


w(r) = - f G 3 a {r - r') /(r') dV, (41) 

P J M 3 

where the Green type function is given by 

11 - 1/2 poo 

G i( r ) = (27r ) 3 / 2 J + c 2 l^| 2 ) A 3 ^i/ 2 (AM)dA. (42) 

Here J 1/2 is the Bessel function of the first kind. 

Let us consider W. Thomson (1848) problem [35] for the fractional elasticity model described 
by equation ( 1371 ) . This problem implies that we should determine the deformation of an infinite 
clastic continuum, when a force is applied to a small region in it. If we consider the deformation 
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at distances |r|, which are larger compared to the size of the region, we can assume that the 
force is applied at a point, i.e. 


/(r) = /o<S(r) = f 0 6(x)6(y)6(z). 


(43) 


Then the displacement held u( r) of fractional elasticity has a simple form of the particular 
solution that is proportional to the Greens function 

u ( r ) = — Ga(r). (44) 

P 

As a result the displacement held for the force that is applied at a point (1431) has the form 


w(r) 


J_Jo_ [°° A sin(A|r|) 
2ir 2 p\r\ J 0 c a X a + c 2 A 2 


(45) 


We can distinguish the following two cases: (1) Weak power-law spatial dispersion with 
cx < 2; (2) Weak power-law spatial dispersion with cx > 2. This is due to the fact that in 
nonlocal elasticity theory usually distinguish the following two cases: (1) Fractional integral 
elasticity (cx < 2); (2) Fractional gradient elasticity (cx > 2). 

For the fractional integral elasticity, the order of the fractional Laplacian is less than the 
order of the term related to Hooke’s law. For the fractional gradient elasticity, the order of the 
fractional Laplacian is greater then the order of the Hooke’s term. 


Fractional integral elasticity model 

The fractional integral elasticity model is described by equation (j37l) with ex < 2 of the form 


c 2 Au(r) - c a ((—A)" /2 w)(r) + - /(r) = 0, (0 < cx < 2). 

P 


(46) 


The order of the fractional Laplacian (—A)"/ 2 is less than the order of the first term related to 
the usual Hooke’s law. Note that the continuum equation (146j) of fractional integral elasticity 
is derived from the lattice equations with weak spatial dispersion in the form (136j) with cx < 2. 
The particular solution of equation (I46j) for the force that is applied at a point (|43l) is the 
displacement held 


u r = 


/o 


2ir 2 p |r| 


A sin(A|r|) 
C 2 A 2 + c a A a 


dX (cx < 2). 


(47) 


Using Section 2.3.1 in the book {43], we can obtain the asymptotic behavior for (1471) for 


where 


i ^ Co(a) 


C'fe(cx) 


C'o(cx) = 


k = 


fo 


j | | (2—o;)(Ai+l)H-l 


Ckiot) = 


2i t 2 p c Q 

k 
2 


foC■ 


2n 2 pc k + 1 J o 


( r ->■ 00 ), 

(48) 


(49) 

2 “)’ 

1 sin(z) dz. 

(50) 
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As a result, the displacement field for the force that is applied at a point in the continuum 
with this type of fractional weak spatial dispersion is given by ■u(r) « Co(a)/|r| 3 ^“, where 

0 < a < 2, on the long distance |r| (§> 1. The asymptotic behavior |r| —> 0 for the fractional 
integral elasticity does not depend on the parameter a. 


Fractional gradient elasticity model 

The fractional gradient elasticity model is described by the equation 


c 2 Aw(r) - c Q ((—A) Q/2 w)(r) + - /(r) = 0, (a > 2). 

P 


(51) 


This may be derived from the lattice model with the fractional weak spatial dispersion in the 
form Q5gp with a > 2. The order of the fractional Laplacian (—A)"/ 2 is greater than the order 
of the first term related to the Hooke’s law. If a = 4 equation (T5T]) become the equation (139]) . 
Therefore the case 3 < a < 5 can be considered as close as possible (a ~ 4) to the usual 
gradient elasticity ()39j) . The continuum equation (l5Tf of fractional gradient elasticity is derived 
from the lattice equations with weak spatial dispersion in the form (j36|) with a > 2. 

The particular solution of equation (13111 for the force that is applied at a point (14511 is the 
displacement field 


u( r) 


fo 

27 r 2 p |r 



A sin(A|r|) ^ 
c 2 A 2 + Cq,A“ 


(a > 2 ). 


(52) 


The asymptotic behavior |r| —» oo for the fractional gradient elasticity does not depend on the 
parameter a. The asymptotic behavior of the displacement field u(r) for |r| —> 0 is given by 


fo T((3 — a)/2) 

2 “ 7T 2 y/TT pc a F(a/2) 


r 


a -3 

5 


(2 < a < 3), 


(53) 


u r 


fo 


0 1—3 /a 3/a • / Q / \ ’ 

lnapc 2 Ca sm(37r/a) 


(a > 3). 


(54) 


Note that the asymptotic behavior for 2 < a < 3 does not depend on c 2 . The displacement 
field w(r) of the short distance is determined only by term with (—A)"/ 2 (a > 2) which can 
thus be considered as a fractional non-locality of the gradient type. 


7 Conclusion 

We demonstrate the close relation between the discrete microstructure of a lattice with weak 
spatial dispersion of a power-law type and the fractional integral and gradient elasticity models. 
We prove that fractional elasticity models can be directly derived from the lattice models with 
power-law spatial dispersion. It has been shown that a characteristic feature of the behavior 
of a fractional non-local continuum is the spatial power-tails of non-integer orders. We assume 
that fractional elastic media (and plasma-like media with power-law spatial dispersion [36] ) 
should demonstrate a universal behavior in space by analogy with the universal behavior of 
low-loss dielectrics in time [44 , 4511461, 47; [48]. Note that the universal behavior in time can be 
connected with the interaction between the particles and the environment [45] . 
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